Re-entrant localization of single particle transport in disordered Andreev wires 
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We study effects of disorder on the low energy single particle transport in a normal wire surrounded 
by a superconductor. We show that the heat conductance includes the Andreev diffusion decreasing 
with increase in the mean free path £ and the diffusive drift produced by a small particle-hole 
asymmetry, which increases with increasing £. The conductance thus has a minimum as a function 
of £ which leads to a peculiar re-entrant localization as a function of the mean free path. 
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Transport phenomena in mesoscopic wires with dimen- 
sions much less than the dephasing length have been 
studied for several decades; their mechanisms are nowa- 
days well understood (see [lly,y| for review). Depending 
on the ratio between the wire length d and the mean free 
path of electrons £ one can separate different transport 
regimes. The ballistic regime holds for small wire lengths 
d ^ £, here the conductance is given by the Sharvin ex- 
pression, G = {e^ /ttK)N , where N ~ (kpa)^ is the num- 
ber of transverse modes, a is the wire radius, a d, 
and kp is the Fermi wave vector. For a shorter £ <C d, 
transport is diffusion controlled with the ohmic behav- 
ior G ^ {e'^ /Trh)N£/d, neglecting the weak-localization 
interference between scattered electronic waves. With a 
further decrease in the ratio £/d, the ohmic dependence 
breaks down due to the localization effects: the conduc- 
tance decays exponentially 1\ when £/d<N-\ 

This textbook picture holds if the transverse confine- 
ment of the electronic states inside the wire is caused 
by an insulating gap in the material surrounding the 
wire, which results in elastic scattering of electrons at 
the wire walls with large momentum transfer (normal 
reflections). In the present Letter we consider another 
realization of a normal-metal wire conductor where the 
electronic states are confined by a surrounding super- 
conducting material. The superconducting gap A out- 
side the normal wire suppresses to zero the density of 
states (DOS) of single-particle excitations for energies 
e < A, thus localizing them in the transverse direction 
within the wire. These states are essentially determined 
by the particle-hole Andreev reflections with low mo- 
mentum transfer at the superconducting/normal-metal 
(SN) boundaries. If the normal reflection processes at 
the SN boundaries can be ignored, we refer to such a 
normal-metal conducting region inside a superconduct- 
ing environment as to an "Andreev wire" . An Andreev 
wire can be connected through bulk normal-metal leads 
to an external measuring circuit. Note that our definition 
of Andreev wire differs from that used, e.g., in Ref. 5] 
where it was applied for a normal conductor in an insu- 



lating environment, connected to superconducting leads. 
A simple way to create Andreev wires is to introduce vor- 
tex lines in a type-II superconductor or to drive a type-I 
superconductor in an intermediate filamentary state by 
applying a magnetic field. Andreev wires can be manu- 
factured artificially in the form of normal channels in a 
superconducting matrix, using modern nano-fabrication 
techniques employed also for producing a wider class of 
hybrid SN structures such as Andreev interferometers 
and billiards 0- Experimentally, the main distinction 
between Andreev wires and usual conductors is that the 
measurements of thermal conductance are more appro- 
priate to probe the single electron transport in Andreev 
wires because the single-particle part of the charge trans- 
port is short-circuited by the supercurrent. 

As shown in Ref.j^ 0| in the ballistic limit £ ^ d, 
Andreev processes suppress the single electron transport 
for all quasiparticle trajectories except for those which 
have momenta almost parallel to the wire thus avoid- 
ing Andreev reflection at the walls. The particles con- 
fined due to Andreev reflections do also participate in 
the transport but through a slow drift along the trans- 
verse modes (Andreev states) with the group velocity 
Vg = h,~^dtk^/dkz ~ e/pF much smaller than vp- This 
Landauer-type drift contribution is the lower limit of con- 
ductance reached as the contribution of freely traversing 
trajectories decreases with increasing d. In total, the elec- 
tronic heat conductance due to these two mechanisms 
is much smaller than what could be derived from the 
Wiedemann-Franz law using the Sharvin conductance of 
a normal conductor: The effective number of conducting 
modes N^s = hn/T is much smaller than (kpa)'^. The 
conclusion of the suppression of the single electron trans- 
port in clean systems is in a good agreement with exper- 
iments on the heat conductivity in type-I and type-II su- 
perconductors in the direction of magnetic field 0, 0| . 

In the present paper we investigate how the low en- 
ergy transport with e A in an Andreev wire of a ra- 
dius a » ^ is affected by a weak disorder introduced by 
impurity scattering. We consider clean wires £ a and 
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neglect inelastic processes assuming ^ I. We start our 
analysis from a qualitative physical picture elucidating 
the main results of our work. Let us consider a quasi- 
particle propagating within the wire along a trajectory 
that bounces from the normal/superconducting walls at 
both its ends. Neglecting the slow drift, the distributions 
of particles and holes are equal at the wall due to the 
Andreev reflection. Without disorder induced scatter- 
ing, the distributions remain equal throughout the wire, 
thus the single-particle transport associated with these 
trajectories vanishes. With disorder, the distributions of 
particles and holes deviate from each other by an amount 
proportional to the probability of scattering, a/-£, accu- 
mulated on their way in between the two walls. In the 
presence of a temperature difference at the ends of the 
wire, the driving force on the trajectory is proportional 
to (a/d)(Ti —T-i) thus the thermal conductance becomes 
K.A = {T/h)NA, where the effective number of modes is 



(1) 



The counter-intuitive behavior of the single-particle con- 
ductance K.A which increases with decreasing i was first 
predicted by Andreev [l^ (see also 0). The coef- 
ficient A in iQ appears to be a slow function of t. 
A ~ ln(f/a) [iJl- For such "Andreev diffusion", disorder 
with a £ d stimulates the single-particle transport 
as compared to that in the ballistic limit: it opens new 
single-particle conducting modes blocked by Andreev re- 
flections in the ballistic limit. This differs from the dis- 
order effects in normal- metal/insulator/superconductor 
systems where disorder opens two-particle tunneling pro- 
cesses for electrical conductance, see j|3| for review. The 
conductance ka reaches its maximum when the mean free 
path decreases down to € ~ a; it further transforms into 
kd ~ {T/h)N£i for a dirty wire £ <^ a where [l^ 



Nd vpa^D/d - [kpafl/d 



(2) 



Here vp is the normal-state DOS at the Fermi level, and 
D = vpt/i is the diffusion coefficient. For a large mean 
free path £ ^ d the conductance ka transforms into the 
ballistic expression ^] k (x d~'^. 

We find that disorder also modifies the single-particle 
transport due to the drift along Andreev states with 
a group velocity Vg, which results from a small non- 
quasiclassical particle-hole asymmetry. The character- 
istic mean free path for the drift appears to be i^s = VgT 
which is considerably shorter than the usual mean free 
path £. The thermal conductance kl — {T/K)Nl as- 
sociated with the disorder-modified drift is found to be 
proportional to the mean free path for IcG ^ d 

Nl ~ vpa^D^s/d ^ {kpaf {vg/vpf {£/d) (3) 

where Deff = Vn£cS = vf,T is effective diffusion coefficient 
much smaller than D that appears in Eq. JS]). This drift 
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FIG. 1: Number of conducting modes (upper curve) normal- 
ized to that in a Sharvin contact, u = N/ (kpa)'^, as a function 
of the mean free path £ interpolated according to Eqs. 
@ for the wire length d ^ aivp/vg). The grey area at the 
bottom is the localization region, TV ~ 1 corresponding to 



saturates at the ballistic Landauer-type expression j^] 

Nl - {kFa)^Vg/vF) (4) 

for very long £ when ^off ^ d. 

The total heat conductance k = {T/h){NA + Nl) in- 
cludes the Andreev diffusion decreasing as £~^ and the 
diffusive drift that increases with increasing £. Equations 
- I0J are illustrated in Fig. ^ as functions of £. If 
d ^ a{vF /vg), the number of modes has a minimum 



{kpaY {a/d) (vg/vp) 



at £min ^ a,{vF/vg). According to the Criterion of Ref. 
a decrease in the number of modes down to N^ia ^ 1 may 
lead to localization of the transport. Thus, varying the 
mean free path around £^in we can obtain a peculiar re- 
entrant localization: for a long enough Andreev wire the 
quasiparticles may become localized not only in a dirty 
limit £ <^ a but also for a quite long £ ^ a in. such a 
way that the conduction opens again through either the 
quasiparticle drift for longer £ or through the Andreev 
diffusion for shorter £. This occurs for a wire length d > 
dc where dc ~ a(kFa)'^{vg/vF)- 

Model. - To develop a more quantitative theoretical de- 
scription we use a quasiclassical approach modified tak- 
ing account of the trajectory drift along the Andreev 
wire. The excitation spectrum is |l5| 



(5) 



Here pz — pf cos 6 and 2xc ^ a is the length of the 
projection of the trajectory section between two walls 
onto the plane perpendicular to the wire axis (the z axis), 
see Fig. 13 Due to a pz dependence of the energy, particles 
perform a slow drift along z with a group velocity 
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FIG. 2; Trajectories for particles (1), (3) and holes (3), (4). 

The group velocity can also be obtained |0| by consid- 
ering particle and hole trajectories, Fig. [3 The momen- 
tum along a; is = ± rne/p_L for particles and holes, 
where p_L = \/v^f ~ Pz — PpsvaO. These trajectories 
have the direction angles tan6'± = p^/pz with respect 
to the z axis such that sin (6*+ — 9-) = (c/cf) cot 0. The 
velocity along the trajectories are 



v± 



±\2 



pI=vf± e/pF ■ 



(7) 



Consider a particle trajectory (1) in Fig. |21 with the co- 
ordinates X = s+sin0+, z — zq + s_|_cos0+, where s+ 
is the distance along the trajectory. The particle is An- 
dreev reflected at the wall and returns back along the hole 
trajectory (2) with the coordinates x = s_sin0_, z = 
z'q + s_ cos 0_ . Both s± are measured from the wire axis. 
The intersection with the wall has the coordinates x = Xc 
and z = zq + Xc cot 0+ = z'q + Xc cot 9- . The drift veloc- 
ity defined as Vg — vf sin9{z'Q — zq)/2xc coincides with 
Eq. (|SJ). This approach is valid for wide wires a S> ^ 
where the drift Zc is much larger than the width kp^ 
of the wave packet. Since Zc ~ a{vg/vF), this requires 
{vg/vF){kFa) 3> 1 which is equivalent to the condition 
e ^ Co, where eo is the lowest level energy in Eq. |5j. 

Kinetic equations - For e ;3> eq the DOS in the normal 
region coincides with the DOS vf in the normal state. 
The Boltzmann kinetic equation is 



dn 1 

V— = [n 
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(8) 



Here (• ■ ■ ) denotes averaging over the momentum direc- 
tions. For particle and hole distributions 

/+(e=P) = np^e, /-(e,p) = 1 - ri-p -e , 
respectively, the Boltzmann equation takes the form 



df± 



1 



(/± " (/±)) 



(9) 



In both the upper-sign and lower-sign equations, the dis- 
tance s is measured in the direction of -t-p. 



Since all particles are Andreev reflected as holes, the 
single-particle current through the wire side walls van- 
ishes. Using mv± sin^-t = p^ we put at the walls 

w+ sin^^+/+ = ti_ sin^?_/_ . (10) 

The kinetic equation on the trajectory (1) gives 

= .f+(0)e-^+/^+ 

Jo 

The function /+(0) is taken aX, x = z — zq. To get 
the corresponding expression on trajectory (2) one sub- 
stitutes /_|_, s_|_, and zq with s_, and z'q, re- 
spectively. Putting s_|_ = Xc/ sm9j^ and s_ — xd sm9_ 
for and respectively, we insert the result into the 
boundary condition Eq. 110|l . The distributions at the 
trajectories (3) and (4) for s < arc found by replacing 
Zq z'q. The boundary condition for them is applied at 
s+ = — Xc/sin^-i- for and at s^ = —xd sin0_ for /_. 
Assuming that (/(x, z)) depends only on z, 



(/W) = (/(O))+scos0 



dz 



the next step is to expand the boundary conditions in 
small Vg/vF and Sc/( for angles 9 3> a/£. Neglecting 
small terms df2/dz we obtain at the wire axis x = 



h = 



Vf cos 9 



fi - cos ( 



2e dz 



(/i - (/i)) 



Vf dz 

We denote Sc — xd sin0 and introduce 

/! = -(/++/-), /2 = -(/+-/-) 



(12) 



(13) 



(14) 



in accordance with the deflnitions used in the theory 
of superconductivity [l^ . The functions /i and /2 are 
nearly constant along the trajectory within the wire. 
Analysis of Eq. (|ll|l shows that for angles 9 ^ ajl, 
the distinction between the usual and the Andreev diffu- 
sion disappears, and /2 — —£cos9{d{fi) /dz) while the 
counterpart of the first term in Eq. H12|) proportional to 
{vg/vF)fi decreases exponentially. 

The first term in Eq. H12|) describes the drift along the 
Andreev states with the velocity Vg. The second term 
is the Andreev diffusion 12]. Equation H13() introduces 
an effective mean free path £cs = {vg/vF)t = VgT much 
shorter than the usual t. In the ballistic limit of very 
long £ such that (.^s 3> d, the distribution /i is constant 
along the wire and 



/2 



VF cos9 



h 



(15) 
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In the most practical limit when i^g <^ d, 

d{fi) 



than d. The total conduction reaches its minimum 



/l = (/l) - VgT- 



dz 



(16) 



and 



/2 = - 



-4 



d{h) 
dz 



vf cos ( 



(17) 

The first term in brackets describes the disorder-modified 
drift. Its relative magnitude with respect to the An- 
dreev diffusion (the second term) is of the order of 
{vg/vpY{£/ Sc)^ , i.e., much larger than non-quasiclassical 
corrections of the order {e/EpY to the usual diffusion. 
We neglect those corrections in what follows. 
The energy current has the form |l7i | 



If. 



d^r 



dflp 
An 



+ CXJ 



cos6'e/2de. (18) 



For very long mean free path £eff 3> d the distribution 
is determined by Eq. (|15ll . and we recover the Landauer 
formula derived in Ref. [9j 



d^r 



dClp 
in 



evg fi de 



(19) 



With /i = tanh(e/2ri) for Vg > and /i = tanh(e/2T2) 
for < we arrive at Eq. |0J for the heat conductance. 

Assuming that the group velocity is independent of 9 
we obtain a qualitative behavior described in introduc- 
tion. Indeed, for i^s ^ d, the distribution obeys Eq. 
(|17|1 where the last term does not contribute to the cur- 
rent. The current becomes l£ = k(Ti — T2) with the total 
thermal conductance k = {T / h){N a + N l) given by Eqs. 
(0 and (PJ. This simplified picture has to be modified, 
however, due to a rapid divergence of Vg at small an- 
gles. This leads to a more complicated behavior of the 
drift contribution to the heat conduction as a function of 
temperature and of the mean free path characterized by 
different power laws in different regions of i and T. The 
detailed analysis of the heat conduction will be published 
elsewhere. Here we discuss its behavior for a long wire 
with d ;» a^/Ep/T. For relatively short mean free path 
a <C ^ ^ a^/EpjT^ the angular integral of the first term 
in the brackets in Eq. I|17|l is cut off by the exponential 
decay of the drift term for 9 ^ a/i. Therefore, 

NL-ikpaYiT/Epfie'/a^d) . 

This b ehavior is replaced with Nl ~ {kpaf{T / Ep){i/ d) 
for a^jEpjT £ d and further transforms into 

Nl - {kpaf{T/Ep)L 

where L - ln(^/rf) for d < ^ < diEp/T) and L - 
ln(£;F/r) for e > d{Ep/T). Therefore, the full satu- 
ration at the ^-independent Landauer-type drift conduc- 
tion Eq. H19|) occurs when l^a ^ {T/ Ep)£ becomes larger 



iVmin ~ {kpaf{a/d)^T/Ep 



at ^min ^ a-\/ Ep /T. The re-entrant localization is possi- 
ble for wire lengths longer than dc ~ a{kpa)'^ y^T/ Ep. 

To summarize, we develop a theory of single electron 
transport and re-entrant localization in clean Andreev 
wires. Our results could stimulate experimental research 
of these phenomena in the mixed or intermediate state, 
as well as in hybrid SN structures. 
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